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Abstract 

We present the classical solutions of the two-center MICZ-Kepler and MICZ-Kepler-Stark systems. Then we suggest 
the model of multi-center MICZ-Kepler system on the curved spaces equipped with so(3)-invariant conformal flat metrics. 

Dh- 1 Introduction 

CD ■ 

As is known, the Kepler system plays a special role among the finite-dimensional classical and quantum integrable systems. 
J"""- ' On the one hand, it has a number of physical applications, on the other hand, this system is one of the key components 
of hidden symmetry concept [T]. Various generalizations of this system were studied during last century. Among them are 
^ (— i the Kepler problems on spheres and hyperboloids [5] as well as their higher dimensional counterparts. Further important 
O |. generalization of the Kepler problem was constructed in the works of Zwanziger, and Mcintosh and Cisneros [3]. They 
^ replaced the Coulomb center of the system by the Dirac dyon (electrically charged Dirac monopole [3] ) , and added to the 
initial Coulomb potential the specific centrifugal term s 2 /2/j.r 2 , where s — eg/c is the so-called monopole number (e is the 
i electrical charge of the probe particle and g is the magnetic charge of monopole), /i is a particle mass. As a result they 
G obtained the integrable system which is very similar to the underlaying Kepler one. It is presently known as a MICZ-Kepler 
system. For instance, MICZ-Kepler system inherent, besides the rotational symmetry, the hidden symmetry of the Kepler 
system given by the Runge-Lenz vector. Similar to the case of pure Coulomb system, the classical trajectories of the 
MICZ-Kepler system are second order curves. The quantum mechanical properties of Coulomb and MICZ-systems have 
also much in common. For example, the only difference in their energy spectra is the lift of the low bound of the orbital 
quantum number from to \s\. 

The MICZ-Kepler system describes the motion of the electrically charged scalar particle in the field of the Dirac dyons. 
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The origin of the additional centrifugal term -^pi can be also understood in this context. In the presence of monopole 
magnetic field the angular momentum of the system get additional, spin-like term, J = r x n — sr/r, where tt is dynamical 
momentum containing the vector potential of the monopole magnetic field. The magnetic field of Dirac dyon is B = gr/r 3 . 
' On the other hand, there is linear relation between orbital and magnetic momenta of single particle, Ai — ^^J- Hence, 
the interaction energy of this magnetic momentum with magnetic field B is given by the expression: 
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i.e. coincides with the centrifugal term in the MICZ-Kepler system [5]- Though, the additional part of the orbital 
momentum proportional to s should rather be assigned to the electro-magnetic field than to the particle, this non-correct 
interpretation nevertheless leads to the correct expression. Let us also mention, that MICZ-Kepler system also describes the 
relative motion of two Dirac dyons with electric and magnetic charges (ei,gi) and (62,32) with eig2 — e2<?i = cs. Similar 
to above mentioned case, the centrifugal term s 2 /2/zr 2 could be interpreted as the interaction energy of the induced dipole 
momentum with the dyon electric field plus interaction energy of the induced magnetic moment with magnetic field (here 
u. is the reduced mass). 

The analogy between MICZ-Kepler and Coulomb systems has elegant explanation in terms of four-dimensional space: 
both systems can be obtained from the Hamiltonian of four-dimensional oscillator via the reduction by U(l) group action 
[5J. In the same way, one can construct the five-dimensional analog of the MICZ-Kepler problem, reducing the eight- 
dimensional oscillator by SU{2) group action [7]. In this case, the SU(2) Yang monopole appears in the system instead 
of Dirac monopole. In two-dimensional case one can also construct MICZ-Kepler like system from the oscillator: for this 
purpose the two-dimensional (quantum) oscillator should be reduce by the action of the parity (Z2) group. The emerged 
two-dimensional MICZ-Kepler system is distinguished by the presence of magnetic flux which provide the system with 
spin 1/2 (Z2 monopole) [8j. The reduction procedure mentioned above have very deep and close relation with the Hopf 
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maps (see, e.g. [9]). For the detailed review of the quantum mechanical aspects of these issues the reader may see the Ref. 
[TU] . The MICZ-Kepler systems on the three-dimensional spheres [TT], on two-, three- , and five-dimensional two-sheet 
hyperboloids [12] as well as their generalization to arbitrary dimensions [13] are also distinguished by their close similarity 
to the corresponding Coulomb systems. For example, the only difference between the quantum-mechanical properties of 
Coulomb and the corresponding MICZ-Kepler systems consists in the lift of the range of the total angular momentum, 
which in its turn leads to the degeneracy of the ground state. Indeed, any three-dimensional rotationally invariant system 
(without monopoles) will preserve its main properties if incorporation of the the Dirac monopole in the center will be 
supplied by the following change of potential 

U(r)^U(r) + ^, (2) 

where G(r) defines the conformal flat metric of the configuration space, ds 2 = G(r)(dx) 2 . Particularly, the form of classical 
trajectory of the system, (p — f(r) unaffected by this kind of replacement, whereas the orbital plane will not be orthogonal to 
the orbital momentum. Instead, the angle between them will be appeared which is defined by the expression cos 9 = s/ J[9 . 
The actual quantum-mechanical spectrum of the extended system coincides with that of the initial one, with minor modifi- 
cation of the possible range of orbital and magnetic quantum numbers j,m: j = |s|,|s| + l,..., m = — \s\,— |s| + 1, ...,|s| |14j . 

The problem of particle moving in the field of two Coulomb centers (or two-center Kepler problem) was solved in the 
middle of XIX century by Jacobi. He established the integrability of two-center Kepler system and of its limit when one of 
the forced centers is placed at infinity which yields the homogeneous potential field (we shall refer it as Kepler-Stark system) 
in elliptic and parabolic coordinates respectively. However, the generalization of this picture to case where Coulomb centers 
are replaced by dyons has been lacking until now. In our recent paper with S.Krivonos [15] we proposed the generalization 
of the MICZ-Kepler replacement $2$ which can be used in case of N Dirac monopoles, 

where — egi, with gi be the magnetic charge of the i-th monopole located at the point with coordinates a^, and rj = |r— a|. 
This replacement has the following important features: 



The system (without monopoles) admitting separation of variables in elliptic/parabolic coordinates results in the 
separable two-center MICZ-system ([3|) with the Dirac monopoles placed at the foci of elliptic/parabolic coordinates. 

The system admits the AT = 4 supersymmetric extension at the following choice of potential 



The corresponding supersymmetric system was constructed in earlier paper by Ivanov and Lechtcnfcld 16J. 

On the Euclidean case the generalized MICZ-Kepler system yields the following generalization of the multi-center Kepler 
problem describing the motion of the particle with electric charge e in the field of N static Dirac dyons 




(5) 



where 



A D r = 6 

r(r — nr) 

is the vector-potential of the Dirac monopole, gi{qt) are the its magnetic(electric) charge of i-th dyon, and is the unit 
vector pointed along the Dirac singularity line. 

Since the two-center Kepler/Kepler-Stark) system admit separation of variables in elliptic/parabolic coordinates, the 
expression mentioned above provides us with integrable two-center MICZ-Kepler system. Moreover, when the background 
dyons have the same ratio of the electric and magnetic charges (i.e. obey trivial Dirac-Schwinger-Zwanziger charge quan- 
tization condition), the potential of the system ([5|) belongs to the class (j4]), i.e. the system allows N — 4 supersymmetric 
extension. It is also noteworthy that the class of integrable models constructed in Ref. [17] include such two-center 
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MICZ-Kepler system in flat space as the limiting case. However, even classical trajectories of the suggested two-center 
MICZ-Kepler systems were not properly studied in the mentioned paper. Another disadvantage of the substitution ([3]) is 
that it is ill-defined for the MICZ-Kepler systems on spheres and hyperboloids. Namely, extending in this way the (one- 
and two-center) Kepler system on the sphere and hyperboloid one will obtain the system which a priori is not separable 
in the same coordinate systems as the initial one. Also these systems do not belong to the class allowing the Af = 4 
supersymmetric extensions. 

In this paper we discuss these issues. 

In the Second Section we describe in detail the two-center MICZ-Kepler and MICZ-Kepler-Stark systems and present 
their classical solutions. 

In the Third Section we present the alternative model of the multi-center MICZ-Kepler system on three-dimensional 
sphere and hyperboloid, as well as on any SO(3) invariant spaces. 



2 Integrability of two-center MICZ-Kepler system 

Now let us demonstrate the classical integrability of the two-center MICZ-Kepler system, and construct the solutions. 
Suppose, that background dyons with magnetic(electric) charges 91,2(91,2) are fixed on z— axis at points (0,0, —a) and 
(0, 0, a) respectively. Choosing the appropriate gauge for the vector-potentials in spherical coordinates 

A^=A\ 2 =Q, A^=. 91,2 cos 0i,2, (7) 

we arrive at the following Hamiltonian: 

n, vl , Pe , ip<p ~ sicosfli - s 2 cos6> 2 ) 2 1 f s x s 2 \ 2 eqi eq 2 



2 2r 2 2r 2 sin 2 6> 2 \n r 2 J n r 2 

In virtue of the relevant form of additional quadratic potential term the system allows separation of variables in elliptic 
coordinates given by the formulas: 

t r ± +r 2 n - r 2 

Using the general kinematic relations resulting from the geometry of problem 

r\ sin 0\ — r sin 9, r 2 sin 9 2 = r sin 9, 

cos(9 1 -9 2 ) = 4a2 - rl - r \ (10) 
2rir 2 

r cos 9 + a r cos 9 — a 

COStq = , COS #2 = ; 

ri r 2 
one obtains the following expression for the Hamiltonian ([5]): 

U = 2a*{?-rf) + (1 ~ ^ + Vi0 + W(V)) ' (U) 

where 

P 2 r, + si - 2p l0 s_f p 2 + s 2 - 2v i0 s+ri 

V(0 = Pv ± + 2aeq + t W r? = ^ ~ ± 1 + 2a eq- V . (12) 

£ z — 1 1 — rj z 

Here q± = q\ ± q 2l s± = si ± s 2 . Now choosing for the generating function the expression S = Lpp v + Si(£) + S 2 (7]) — Et 
we arrive at the following Hamilton-Jacobi equations: 



(13) 



(C 2 1) (^) 2 + V(0 - 2a 2 E (e - 1) = n, (l - r, 2 ) (j^f + W M - ^E (l - v 2 ) = -n, 
where n is the separation constant which fixes the values of the corresponding constant of motion 

= -^2 (« 2 - M + V (0) + e ((1 - V 2 )P 2 , Mr?))) (14) 
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Then, integrating Eqs. (|13[) one obtains the generating function of the system 



S = p v <p~Et + J J 2a*E + ^— - (g2 _ 1)2 ^ + y W 2a 2 ^ — ^ (1 - ^ d?7 (15) 

Hence, we proved the exact solvability of classical two-center MICZ-Kepler system. Now, taking derivatives of the 
generating function by constant of motion, dS/dp v , dS/dn and dS/dE and putting them to constants (which we are 
claiming, without loss of generality, to be zero) one obtains the classical solution of the system. Its trajectories are defined 
by the expressions 

dt, f dr, 



(e-i)^ET^f J (i-^)^e-^i 

(p v - s-g)d4 f (p v - s + rj)dr] 



The time evolution of the system is given by the expression 



t = « 2 / , i +« 2 / , dv as) 



Here t is the time and V(£) and W(?7) are defined in Eqs. (fT2|) 
2.1 MICZ-Kepler-Stark(-Zeeman) system 

The system introduced above has an important limiting case. Putting one of the dyons at the infinity one obtains the 
Hamiltonian of the charged particle moving in the filed of one localized dyon with additional parallel homogeneous (constant 
uniform) electric and magnetic fields: 

(p-eAn-eBxr/2) 2 1 /s „ \ 2 eq „ „„„ , , 

H = ^- — n — + -+ eBr + — -eEr, EB 19 

2 2 \r / r 

Due to the presence of constant uniform magnetic and electric fields one could refer this system as MICZ-Kepler-Stark- 
Zeeman system. As in the previous case, here one have the integrable system, i. e. the system which allows separation of 
variables in parabolic coordinates. 

In order to show it we, at first, turn to spherical coordinates in which the vector potential of the homogeneous magnetic 
filed assumed to be in z-direction is 

A r = A e = 0, A v = ^Br 2 sin 2 9 (20) 

thus, the Hamiltonian takes the form 

H = 4 + A + iP * - C ° S ° : H 13 ^ Sln2 9) 2 + \ C 1 + eBz) \^~eE z (21) 
2 2r 2 2r 2 sin 2 6» 2 \r J r K ' 

Due to additional quadratic term arising from the relevant MICZ-Kepler replacement the corresponding Hamilton-Jacobi 
equation become separable in parabolic coordinates given by the following relations: 

£ = r + z, n = r-z (22) 

The Hamiltonian then reads 

n = (4£p 2 + 4 V p 2 + V(0 + W(t])) - \ Pv eB, (23) 



where 



V(0 = {Pv j + 3seBH - eE^ 2 + + 2 + eq, W(t?) = {P<P - ^ - ZseB^ + eEif + + 2eq. (24) 
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Then, taking for the generating function S = <pp v + <Si(£) + $2(11) — (E — ^p v eB)t one arrive at the following Hamilton- Jacobi 
equations: 

(dS x \ 2 , ,„ >x _ , /dS 2 



4 H^y +v(o-2^e = n, M^fJ + M*?) - 2Ei i = - n > ( 25 ) 

Here, us usual, n is the separation constant which fixes the values of specific constant of motion, responsible for the 
separation of variables. For this quantity from Eqs. (f2"3")l and (f2"5)) one obtains 

Ip = j^r (^V(P 2 £ - V%) + »?V(0 - W?)) (26) 
The generating function obtained from the corresponding Hamilton- Jacobi equations is 



5r=^_(£?_ip v B)t+i | ^E-iseB-^^ + ^^+eEi-^edi- 



\JJ 2 E + ZseB-^- S ?-^ 



2£ + 3sei3 - ^ y _ 2 y - - eity - —j-ifd V . (27) 



Hence, we prove the classical integrability of the MICZ-Kepler-Stark-Zeeman system 

From the generating function we immediately get the expressions for the trajectories of system 



/ ^ = I = (28) 



1 Di , 1 f (p v + s)d£ 1 f (p p — s)d?7 



u = -fli + - / ^ ' y s - + - / — f v ' ' (29) 

2 27 t /op , n-Vfjy 27 „ / OF n+WW 



and for the time evolution 



Here f is the time and V(£) and W(r?) are defined in Eqs. ([24 



1=1 * + [ . dV (30) 

2E + J j2E- n+WM 



3 Multi-center MICZ-Kepler system on curved space 

As was mentioned in Introduction, any system (on conformal-flat space) without monopolcs admitting the separation 
of variables in elliptic/parabolic coordinates remains separable in the same coordinate systems under incorpration of 
monopoles supplied by the corresponding modification of potential term given by ([3]). The two-center MICZ-Kepler 
systems considered in previous Section are particular examples of these systems. However, this procedure fails in case of 
the Coulomb systems on spheres and pseudospheres (two-sheet hyperboloids). Due to the existence of hidden symmetries 
given by Runge-Lenz vector the Coulomb system (on Euclidean space) admits the separation of variables in spherical, 
elliptic and parabolic coordinate. The Coulomb system on (pseudo)spheres also has the hidden symmetries given by the 
analog of Runge-Lenz vector. Analogously to Euclidean case, this hidden symmetry is connected to the the separation 
variables in few coordinate systems. These coordinate systems turn to the spherical, elliptic and parabolic ones at the 
planar limit. The connection between spherical and Cartesian coordinates on the sphere is identical with that on the 
Euclidean space, whereas the discrepancy is appeared for the elliptic and parabolic ones (see, e.g., [H] and refs therein). 
Therefore, the statement of the previous Section concerning the separability of variables for "MICZ-extended" systems is 
no more valid for the Coulomb systems on (pseudo)sphere. Moreover, it seems that the substitution §5§ is ill-defined on 
non-Euclidean spaces including (pseudo)spheres. The evidence of it stems out from the following. When the magnetic and 
electric charges of the background dyons obey trivial Dirac-Schwinger-Zwanziger quantization condition 

9i<lj - 9j1i = 0, (31) 
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the potential of the Euclidean multi-center MICZ-Kepler system belongs to the class (01 (up to unessential constant), i.e. 
it admits TV = 4 supersymmctric extension. While the potential of the multi-center Coulomb system on (pseudo)sphere, 



1 1 



U(r) = e<f) q i,.., 9f( =ey'?^(r i ), where cf)(r) = , ri = |r-aj| (32) 



i=l 



does not belong to the class (UJ). Here e = 1 correspond to the sphere, and e = — 1 to the pseudosphere, and the metrics 
looks as follows 

2 _ 4r 2 (rfr) 2 

The (one-center) Coulomb potentials on Euclidean space, sphere and pseudosphere are nothing else, but the so(3) invariant 
Green functions for the Laplasians defined by the corresponding metrics, A</jc = S(r). While the vector potential of the 
Dirac monopole is just a one-form dual to this Green function. Explicitly, 

*dA D = -d4> c , =>• A</j c = S(r), A = *d*d + d*d* (34) 

Taking into account the duality between vectors and one- and two-forms in the three-dimensional spaces, one can write 
down these expressions in the following way 

- V<fc = rotA^, =► A<j) C = 5(r). (35) 

Let us write down the explicit so(3) invariant solutions of this equations on the so(3) invariant space with the metrics 
ds 2 = G(r)(dr) 2 . The (zero-, one and two-)form fields are independent on the metric, hence, one has to choose 

. , (r x dr) A dr . d<h(r) rdr . , 

dA D (r) = i ^ , #(r) - -g^ — . (36) 

Taking into account the conformal flatness of the metrics, we get, from (f3~4"]) 

d(f>c 1 i f df 



dr r 2 ^G^r) J r 2 ^G^r) 

Particularly, for the Euclidean space one obtains <j>c — 1/r, and for the (pseudo)sphcrc (f3"3"| the potential given in (J32J) 
Now, we define the multi-center MICZ-Kepler system on the so(3) invariant space by the expression 



(37) 



(p eA gi! ..., gjv ) 2 e 4> gi ^. a< 
H— 1 o re<p gu ... t q N , (38) 



where A.,, ,. x = 5iA D (r - a,), (f) qi ,..., qN = Ei%^c(n)- 

The corresponding one-center Hamiltonian on the sphere and two-sheet hyperboloid coincides with the MICZ-Kepler 
Hamiltonian constructed within the rule ^ up to unessential constant es 2 /4r 2 . Hence, it is separable in the spherical 
coordinates, and in the modified elliptic and parabolic coordinates considered in Ref. 18] . It seems, that two-center 
MICZ-Kepler systems on spheres and pseudospheres are also separable in these coordinates. 

When the electric and magnetic charges of background dyons obey the condition (|31[) . or, equivalently, q%/gi = en, this 
Hamiltonian can be represented in the form 

n, : P 'A,,. , x i 2 e 2 (0 3l ,...^ 3JV +K/e) 2 k 2 

n = 2G + 2 y (39) 

where grad(</> Slj ... igjv + n/e) = — rotA ffl> ... i9N . This Hamiltonian admits the TV = 4 supersymmetric extension on the Eu- 
clidean space ([16]) and on the sphere ([19]) (provided the unessential constant k 2 /2 is omitted). In our knowledge, the 
supersymmetric extensions of this systems on the generic conformal-flat space are unknown. However, we believe, that 
with the explicit component expressions for the supercharges of the system on sphere [?], we will be able to guess the 
supersymmetric system on any conformal flat case. 
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